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ABSTRACT 

In energetic nonthermal sources such as gamma-ray bursts, AGN or galactic 
jet sources, etc., one expects both relativistic and transrelativistic shocks 
acompanied by violent motions of moderately relativistic plasma. We present 
general considerations indicating that these sites are electron and positron 
accelerators leading to a modified power law spectrum. The electron (or e^) 
energy index is very hard, oc or fiatter up to a comoving frame break energy 
7^,, and becomes steeper above that. In the example of gamma-ray bursts the 
Lorentz factor reaches 7,, ~ 10^ for accelerated by the internal shock ensemble 
on subhydrodynamical time scales. For pairs accelerated on hydrodynamical 
timescales in the external shocks similarly hard spectra are obtained, and the 
break Lorentz factor can be as high as 7^, ^ 10^. Radiation from the nonthermal 
electrons produces photon spectra with shape and characteristic energies in 
qualitative agreement with observed generic gamma-ray burst and blazar 
spectra. The scenario described here provides a plausible way to solve one of the 
crucial problems of nonthermal high energy sources, namely the efficient transfer 
of energy from the proton fiow to an apropriate nonthermal lepton component. 
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1. Introduction 

Shocks, and more generally systems of shocks as well as turbulent flow downstream 
from such shocks may be a common feature of a number of nonthermal gamma-ray sources. 
A recent example of interest is the dissipative relativistic fireball model of gamma-ray 
bursts (GRB), which follows from very general energetic and observational constraints, 
independently of whether the sorces are at cosmological or galactic halo distances. This 
implies a highly relativistic outflow of matter and electro-magnetic flelds lasting on the order 
of seconds where both external and internal shocks are expected (e.g. Rees & Meszaros, 
1992, 1994, Narayan, et.al, 1992, Paczyhski& Xu, 1994), providing a scenario in good 
qualitative agreement with major observational requirements (e.g. Meszaros, et.ai, 1994, 
Meszaros& Rees, 1994). Under the cosmological assumption (which we henceforth assume) 
likely energy sources may be, e.g. the coalescence of a compact binary (Narayan, et.ai, 
1992) or a failed supernova-hke collapse (Woosley, 1993). While the blast wave propagating 
into the external medium is highly relativistic, the reverse shock propagating back into the 
ejecta is likely to be only moderately relativistic. In addition, the irregular nature of the 
primary energy release results in the formation of a complex internal structure of faster 
and slower portions of the flow leading to internal shocks having moderate Lorentz factors 
(~ 1) in the comoving frame of the wind. Hydrodynamic (Waxman & Piran, 1994) or MHD 
(Thompson, 1994) turbulence may be expected in such scenarios and could play a role in 
the flow dynamics. In AGN jets also one expects internal shocks, as well as termination 
shocks, and interestingly, the gamma-ray spectrum of blazars is qualitatively similar to 
that of GRB. Such systems of shocks and turbulent regions provide an environment similar 
to those thought to lead to efficient particle acceleration (e.g. Blandford & Eichler, 1987, 
Jones & Ellison, 1991, for a review). 

In the case of GRB, which we take as a generic example in this paper, the shocks 
are expected to energize an interaction region of spatial scale A ~ ct^ar^ in the wind 
comoving system, where in the case of internal shocks t^ar is the time scale of energy 
release fluctuations (or in the case of the reverse external shock it is the light crossing time 
over the energy deposition region) and F is the mean bulk Lorentz factor of the flow. A 
multiple-shock structure is likely to arise as a result of the energy release fluctutations or 
the reflection and intersection of shocks crossing flnite shells or an inhomogeneous outflow. 
Strong, smooth relativistic MHD fluctuations produced by the irregular flow motions are 
expected to arise and interact with particles accelerated in such shocks. While details are 
uncertain, it is possible to consider some fairly general features of such flow collision regions 
(PGR), encompassing an ensemble of internal shocks and developed turbulent motions, of 
maximum comoving lenght scale A and occuring at a lab frame distance ~ ct^ar^^ from 
the center of the event. Large-amplitude variations of the bulk velocity and the magnetic 
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field can be expected on comoving length scales I — aA, with a < 1 inversely proportional 
to the number of shocks. For the multiple secondary shock generation expected from flow 
collisions we estimate a < 10~^. 

In addition to such larger scale hydrodynamic substructure, we suggest that instabilities 
and nonhnear effects provide some sort of cascading process leading also to a wide spectrum 
of MHD fluctuations on smaller scales. Note that small scale MHD motions may have 
nonrelativistic bulk velocities (unlike the large scale motions containing most of the energy). 
Such violent systems arc favorable sites for nonthermal particle acceleration. The nearest 
nonrelativistic analogue are the corotation interaction regions of the solar wind, which 
have been known for a long time to be an efficient source of low energy cosmic rays (e.g. 
McDonald et al. 1974). Similar systems were suggested by Rees (1987) as sites for energetic 
particle acceleration in AGN jets. 

In the next section we consider the process of nonthermal electron and positron 
spectrum formation in FCRs. While the specific examples and numerical values refer to 
cosmological gamma ray bursts, similar considerations apply to galactic halo burst models, 
and may also be of relevance for other astrophysical (e.g. AGN or galactic jet) sources 
involving relativistic flows and shocks. 

2. Nonthermal lepton spectra 

We consider here a nonequelibrium processes of transformation of the power in the 
baryon bulk motions of the relativistic wind to nonthermal leptons, and the temporal 
evolution of the lepton spectra. Turbulent plasma motions and shocks have been considered 
as generic sources of nonthermal particles since Fermi's pionering work on statistical 
acceleration (e.g. Blandford & Eichlcr, 1987). There are three important time scales in 
our problem. The first one is of the order of the cyclotron gyration time of relativistic 
e"*^, which is characteristic of the fast prcaccclcration process taking place at shock fronts. 
The spectra of leptons accelerated by single shocks may be steeper than what is needed to 
explain directly gamma-ray burst photon spectra. Thus we consider the individual shocks 
as injection agents providing superthermal particles that can be further subject to diffusive 
acceleration through scattering on resonant fluctuations and large scale MHD plasma 
motions. The latter, as shown below, can produce very hard lepton spectra on longer, 
subhydrodynamical and hydrodynamical timescales. 
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2.1. Cyclotron time scale 

Superthcrmal particles can be naturally extracted from the thermal pool by collisionless 
shocks, and this process operates on cyclotron timescales. The microscopic physics of 
particle acceleration in relativistic shocks is very complicated, and several fundamental 
aspects remain unclear, although some important results highlighting the distinctive 
features of such shocks have been obtained, e.g. Hoshino et.al.{1992). These authors find 
that for transverse relativistic shocks pair acceleration to nonthermal energies occurs if the 
upstream flow contains ions carrying most of the energy flux, and they speculate that this 
may extend also to electron-proton plasmas. (Such conditions are typical of GRB, where 
a proton component is expected and pairs are present but do not dominate the energy 
density). Hoshino, et.aL{1992) obtain a downstream nonthermal pair distribution oc 7^^, 
where N is the number of particles within the range d^, and their results indicate that a 
fraction 0.1 — 0.2 of the upstream baryon flow energy goes into magnetosonic waves which 
accelerate the nonthermal pairs. Microscopic simulations of relativistic quasiparallel shock 
are not yet available, but for nonrelativistic shocks in proton-electron plasmas an injection 
fraction ( ~ 10~^ is typical (e.g. Giacalone et.al, 1992). 

2.2. Subhydrodynamical time scale 

A common attribute of any Fermi-type acceleration process is the isotropization of 
the fast lepton distribution due to the scattering by magnetic field fluctuations. Prom 
momentum conservation one infers that to provide efficient particle scattering the magnetic 
field fluctuation must have an energy density comparable to that of the fast particles in 
the FOR rest frame. Resonant scattering of superthermal leptons will also be accompanied 
by stochastic acceleration with a typical time scale r** cx {c/uphY x (A/c), where Uph is 
the phase velocity of the waves resonating with the scattered particle. For the mean free 
path A (7) of a relativistic charged lepton within such a system of superthermal particles, 
the fluctuation spectrum can be approximated as being continuous on scales much larger 
than the electron gyroradius, due to strong dissipation. For a broad range of magnetic field 
fluctuations with spectral energy density W{k) oc /c""'*, one expects A oc Z x {rg/l^'^'^^ from 
standard quasilinear theory (e.g. Blandford & Eichler, 1987), where Vg — 1.6 x 10^ 7 cm 
is electron gyroradius. The use of quasilinear theory for the effects of fluctuations on scales 
< / is justified because of the relatively small amplitudes of the resonant fluctuations (see 
also Hoshino et.ai, 1992) in the sub-hydro dynamic MHD regime. Thus here we calculate 
the temporal evolution of the spectrum using a standard Fokker-Planck treatment. 

This process leads to a nonthermal electron (or e^) spectrum of the form A^(7) oc (n'y^~^, 
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which develops on a sub-hydrodynamic time scale of the order of a few r^*, which for 
Uph ~ c is 

~ {r./iy-'' X l/c . (1) 

Typically {rg/Vf"'^ <^ 1 for the energies of interest (except for the case 11 — 2). Prom energy 
conseravtion, this energy spectrum extends up to 



7* 



7i mp/me e C , (2) 



where 7j is the initial proton (lepton) Lorentz factor, e < 1 is the portion of the total 
upstream power in baryons pumped into turbulent fluctuations, C, is the lepton injection 
fraction. The difference between the steeper 7"^ spectrum injected near the shock and the 
harder 7^^^ spectrum produced by the MHD ffuctuations arises because of the different 
spatial extent of the acceleration regions. In the former the acceleration and escape time 
scales are comparable, ta ~ n/v?, where u is flow velocity and /t ~ t>A is diffusion coefficient 
with i> ~ c the particle velocity, while tesc ~ 5"^ / k where 5 ~ {v/u)\ ^ n/u is the width 
of particle acceleration region near the shock, so tesc ~ u'^S'^/k,'^ ~ 1. For scattering by 
MHD fluctuations between the shocks, however, the acceleration time is much shorter than 
the escape time: ta ~ {rg/lY~^{l/c) ~ A/c, while tesc ~ A^/k ~ (A/A)^(A/c) > (A/c). 

Since typical MHD turbulent spectra have indices ji with 1.5 < < 2, the resulting 
particle spectra are oc 7"^ or flatter. We conclude that a substantial portion of the power 
in turbulence and some portion ~ e of the upstream baryon power in the flow is transferred 

to the charged leptons at energies near 7*. An estimate of e and would depend on 
uncertain details such as the ratio of the matter to antimatter content in the flow, wave 
modes supported by the plamsa, etc. In the absence of other information, we will adopt 
here C, ~ 10"'^ and e ~ 0.1, compatible with the numerical results of Giacalone, ei(.a/.(1992) 
and Hoshino, ei.a/.(1992). Taking as an example = 1.5, an estimate of the break energy 
for the above typical parameters of the model gives 7* ~ 3 x 10^. 

Synchrotron losses of relativistic pairs have a time scale r^yn ~ 5 x 10^-B^^7^^ s, if B is 
measured in G. Thus for the formation of a hard branch of the pair spectrum at least up 
to the Lorentz factor 7* due to resonant acceleration the scale I must satisfy the condition 
I < 10(3m+i2)/(m-i) ^(m-3)/(m-i) cm. This imples, for S = 10^ G and 7* ~ 1 - 2 x 10^ 

the scale / < 10^ cm if turbulent fluctuations have index ji —2, and I < 10^^ cm\l ji — 1.5. 

On the subhydrodynamical time scale the distribution of nonthermal pairs for 7 7* 
will be highly intermittent, with nonthermal pairs at these energes concentrated in the 
vicinity of the shocks, since they lose their energy before being mixed within the FCR. The 
synchrotron photon spectrum of the system beyond the break might be dominated by the 
brightest spots from some particular shock or the superposition of contributions from a few 
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shocks. The resulting spectral shape just beyond the break may be rather complicated, 
reflecting with some modifications the injection spectrum. 



In addition to the above effects, one can also expect acceleration from processes occuring 
on the longer hydrodynamical comoving time scales of the order of l/c. The electric fields 
induced by turbulent motions of plasmas carrying magnetic fields on different scales lead 
to statistical energy gains of the superthermal charged particles. For nonrelativistic MHD 
turbulence the particle energy change over a turbulent correlation length (or correlation 
time) is small, because the induced electric field is smaller then the entrained magnetic field. 
However, the distinctive feature of statistical acceleration in the relativistic MHD turbulence 
and shocks on larger scales expected in the FCR, is the possibility of a substantial particle 
energy change over one correlation scale, because the induced electric fields are no longer 
small. In this case a Fokker-Planck approach cannot be used. Instead, we argue here that 
it is possible to calculate the energy spectra of nonthermal particles within FCRs using an 
integro-differential equation which is a generalization of the Fokker-Planck approach (for 
details see the review by Bykov & Toptygin, 1993, hereafter BT). 

Consider charged test particles interacting with a wide spectrum of MHD fields and an 
internal shock ensemble produced by the colliding fiows within generalized FCRs. In the 
wind comoving frame, we can assume the fluctuations on all scales up to ~ A (including 
the internal shock ensemble) to be nearly isotropic (for the latter, it is enough if they are 
forward-backward symmetric). The small mean free path A of the superthermal particles 
leads to their isotropy in the frame of the local bulk velocity fluctuations. The assumed 
statistical isotropy of the bulk velocity fluctuations in the comoving frame of the wind 
results then in a nearly isotropical particle distribution, after averaging over the ensemble 
of internal shocks and accompanying motions on scales ~ I. 

To calculate the spectrum of nonthermal leptons accelerated by the ensemble of internal 
shocks and large-scale plasma motions in the FCR (averaged over the statistical ensemble of 
large scale motions) we use a kinetic equation for the nearly-isotropic distribution function 
N — 7^F, which takes into account the non-Fokker-Planck behavior of the system (see 



2.3. Hydrodynamical time scales 



Bykov 1991 and BT), 



dF{T,U) 



g.(0+rood6 A(e - 6) AF(r,ei,t) 



dt 




(3) 
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Here ^ = ln(7/7j), 7, is the Lorentz factor of the injected particles, Qi{C) oc (cnl'^ is the 
rate of nonthermal particle injection, n is the lepton number density in the FCR comoving 
frame. The kernels of the integral equation Eq.(3) determining the spatial and momentum 
diffusion are expressed through correlation functions describing the statistical properties of 
the large scale MHD turbulence and shock ensemble. Following the renormalization method, 
the Fourier transforms of the kernels Df{s) and -D|^(s) are solutions of a transcendental 
algebraic system of equations of the form — $1,2 (-Df, s). Here s is a variable which 
is Fourier conjugate of ^. Equation (3) and the renormalization equations are valid only for 
particles with sufficiently small mean free paths A(7) ^ A. 

The crucial point is that the solution of equation (3) has a universal behavior, only 
weakly dependent on the complicated details of the turbulent system. The stationary 
solution to Eq.(3) with a monoenergetic injection rate Qi has an asymptotical behavior of 
a power-law form, N oc Qa'" (Bykov, 1991), where a = -0.5 + [2.25 + eDl(0)D^^(0)]°■^ 
and we took 6 ~ (Z/A)^. For conditions typical of developed turbulence, the ratio of 
the rate of the scattering to the acceleration rate is Di(0)D^^(0) < 1 (see BT), and for 
< 1 we obtain c ~ 1. This hard 7"^ spectral behavior arises because the acceleration 
time ~ //c ~ aA/c is much shorter than the escape time at the relevant energies, 
Tesc ~ A^/k ~ A^/(/c) ~ A/(ac). The power needed to produce such a spectrum of 
nonthermal particles increases oc jmax, so it is important to understand its temporal 
evolution. 

In the test particle limit, where the backreaction of the accelerated leptons on the 
energy-containing bulk motions is negligible, we have A^(7,t) oc C^7~^ for 7 < 7*(t), where 
7*(t) = 7iexp(t/r„^) and 

oc l/c ~ a(A/c) , (4) 

is the typical hydrodynamical acceleration timescale (see e.g. BT), with ~ few, a < 1 
and A the comoving width of the region energized by shocks. From the energy balance 
equation, when the value 7*(t) ~ ji mp/nie e (^"^ is reached the growth must saturate, and 
the resulting spectrum consists of two branches. One is the hard spectrum A'"(7) ~ ^^7"^, 
for 7 < 7yt, where 

7^ ~ 7j mp/nie € . (5) 

For the typical values of our problem 7^ ~ 1 and ( ~ 10~^ so 7* ~ 10^ (but it could be even 
larger since e ~ 1 for large scale plasma motions). 
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3. Electron Energization Efficiency 

We outline here the apphcation of the above acceleration scenario to the dissipative 
fireball model of gamma-ray bursts. It has been argued in the introduction that one can 
expect violent flow collision regions (FCRs) to form in the dissipative portion of the fireball 
evolution. We do not go here into a detailed discussion of the radiation physics nor do we 
attempt to model GRB in detail, concentrating rather on broadly generic examples using 
typical values of the relevant physical quantities. 

In the fireball wind models, FCRs might occur around radii r ~ ct^ar^'^ ~ 10^^ — 10^^ 
cm with bulk Lorentz factors F ~ 10^ and mean comoving field strengths B ~ 10^54 G (e.g. 
Rees & Meszaros, 1994). We assume the acceleration to occur beyond the region where 
significant pair formation is expected, i.e. outside the photosphere, and take as numerical 

examples a lepton injection fraction ( ~ 10~^, turbulence energy fraction e ~ 10~^, and 
initial injection Lorentz factor 7j ~ 1. We assume also a broad spectrum of MHD or 
whistler type fluctuations with index fi = 1.5. The characteristic time scale of the fast 
injection process of leptons is a few microseconds in the FOR comoving frame. Then, on 
a subhydrodynamical time scale r^* ~ 10c~^(/rg)^/^ ~ 10c~^{actyarTrgy^'^ a hard lepton 
spectrum is established, with N{'-f) oc C^7~'''^ for 7 < 7* ~ [7, nip/nie e C"^]^''^- The 
comoving frame time r^* ^ few ms is enough to transfer a fraction ~ e of the baryonic power 
to the accelerated leptons, with Lorentz factors 7 ~ 7* ~ 2 x 10^. For steeper turbulent 
fluctuation spectra (// ^ 1.5) the timescale is somewhat longer but still compatible with 
observational requirements. As they are accelerated, the leptons radiate a synchrotron 
spectrum of the form uFi, ~ u^''^'^ up to a break near 0.1 MeV for the parameters used, and 
cx u'' above that with k between and -1 if rather weak shocks dominate in the FCR. 

Nonthermal lepton acceleration on the longer, hydro dynamical timescales might be 
important in GRB external shocks around radii r ~ 10^^ cm, e.g. Meszaros, Rees & 
Papathanassiou, 1994. Typical magnetic flelds in the reverse shock could be ~ 10 G, 
while the hydrodynamic time scale is ~ aA/c ~ tens of seconds in the FCR comoving 
frame (and F"^ times shorter in the lab frame). This is enough to form a pair spectrum 
N{'y) oc C^7~^ up to Lorentz factors 7^ ~ 2 x 10^, for these parameters, and a synchrotron 
spectrum with uF^, ~ v and a peak energy near the 0.5 MeV range. 

Recent GRO observations of blazars (McNaron-Brown et.ai, 1995) show clear evidence 
for broken power-law spectra peaked in the MeV range with a shape similar to gamma-ray 
burst spectra (e.g. Greiner, et.al., 1995). From the similarity of the acceleration scenarios 
expected in both, involving relativistic shocks and turbulence, one may speculate on the 
possible applicability of the above processes to explain blazar spectra. 
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In principle, an injected fraction C < 1 of protons may also be accelerated by 
the same mechanisms. The maximum proton Lorentz factors would in this case be 
7p,* ~ (e7j/C)^^''^~^^ or 7p,* ~ (^7i/C)i which for e ~ 10~^, ( ~ 10~^ could be as high as 
10^. However, the fraction of postshock proton energy going into such a fiat, nonthermal 
relativistic proton component is at most e, comparable with the fraction of energy going 
into the nonthermal flat lepton spectrum. 

The efficiency of the transfer of energy from the proton to the lepton component in 
these models is high, typically of order e^e ~ (C/7j)(^e/'"^p)7f^^ ^ e for subhydrodynamic, 
and Cpe ~ (C/7i)(?Tie/'nip)7e ^ e for hydrodynamic acceleration, where 7e can go up to 7^, 
or 7^, and e ^ 1 is the fraction of upstream proton energy converted into turbulence in 
the semirelativistic wind and reverse blast wave shocks, ji is initial Lorentz factor and C 
is lepton injection fraction. This high efficiency is due to the very hard lepton spectra 
achieved with tacc ^ tesc, ^il) 1"^ or flatter, which puts most of their energy near the 
upper break value 7* or 7*. This is a signiflcant fraction of the equipartition value between 
the accelerated Icptons and the bulk of the shocked protons. Similar shocks and turbulent 
regions are likely to be present in AGN or galactic jets. If these are electron-proton jets, 
as opposed to electron-positron jets, most of the energy is in the protons (as for the GRB 
case), and a signiflcant fraction of it should be channeled into the electrons. Since the 
leptons, due to their smaller mass, are responsible for most of the radiation, this mechanism 
of proton-electron energy sharing fulflUs a major prerequisite for a high radiative efficiency 
in GRB and other nonthermal gamma-ray sources. 
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